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ABSTRACT 

We present an analysis of stellar distributions for the young stellar clusters GGD 12-15, 
IRAS 20050+2720, and NGC 7129, which range in far-IR luminosity from 227 to 5.68 xlO 3 L Q and 
are all still associated with their natal molecular clouds. The data used for this analysis includes ncar- 
IR data obtained with FLAMINGOS on the MMT Telescope and newly obtained wide-field 850 //m 
emission maps from SCUBA on the JCMT. Cluster size and azimuthal asymmetry are measured via 
azimuthal and radial averaging methods respectively. To quantify the deviation of the distribution of 
stars from circular symmetry, we define an azimuthal asymmetry parameter and we investigate the 
statistical properties of this parameter through Monte Carlo simulations. The distribution of young 
stars is compared to the morphology of the molecular gas using stellar surface density maps and the 
850 fim maps. We find that two of the clusters are not azimuthally symmetric and show a high degree 
of structure. The GGD 12-15 cluster is elongated, and is aligned with newly detected filamentary 
struc ture at 850 um. I RAS 20050+2720 is composed of a chain of three subclusters, in agreement 
with iChen et alJ (j!997f) . although our results show that two of the subclusters appear to overlap. 
Significant 850 /im emission is detected toward two of the subclusters, but is not detected toward the 
central subcluster, suggesting that the dense gas may already be cleared there. In contrast to these 
two highly embedded subclusters, we find an anti-correlation of the stars and dust in NGC 7129, 
indicating that much of the parental gas and dust has been dispersed. The NGC 7129 cluster exhibits 
a higher degree of azimuthal symmetry, a lower stellar surface density, and a larger size than the other 
two clusters, suggesting that the cluster may be dynamically expanding following the recent dispersal 
of natal molecular gas. These analyses are further evidence that embedded, forming clusters are often 
not spherically symmetric structures, but can be elongated and clumpy, and that these morphologies 
may reflect the initial structure of the dense molecular gas. Furthermore, this work suggests that 
gas expulsion by stellar feedback results in significant dynamical evolution within the first 3 Myr of 
cluster evolution. We estimate peak stellar volume densities and discuss the impact of these densities 
on the evolution of circumstellar disks and protostellar envelopes. 
Subject headings: pre-main sequence — stars: formation — infrared:stars 



1. INTRODUCTION 

It is now commonly accepted t hat most stars form 
in clusters of 100 or more stars l)Lada fc Ladal 120031: 
iPorras et alll200.1 lCaroenterl l2000). In addition, young 
OB stars nearly always a ppear in groups and clus- 
ters of lower mass stars l)Testi. Palla. fc Nattal 119991 
iMegeath et al.l 12002]) . Clearly, understanding the pro- 
cess of forming stars in clusters is of considerable impor- 
tance. Although we have an increasingly detailed picture 
of relatively isolated star formation in nearby dark clouds 
such as those in the Taurus complex, our understanding 
of formation in clusters has been hindered by greater ob- 
servational and theoretical challenges. 

Observationally, the challenges of studying embedded 
clusters are their distance, their spatial density, and their 
association with high column density molecular clouds. 
While the Taurus complex provides a number of exam- 
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pies of isolated star formation at 1 60 pc, most clusters 
are at distances of 300 pc or greater i)Lada fc Ladall200'3T: 
IPorras et~a l. 2003), and thus require greater angular res- 
olution and sensitivity to study. The stars in these re- 
gions may be closely spaced, again requiring high angular 
resolution to resolve the individual stars. Finally, many 
clusters are deeply embedded, and are thus obscured to 
observation at visible wavelengths. For this reason, the 
importance of young stellar clusters was not realized until 
infrared cam eras with detector arrays became available 
ijLadal fl99^ . With the latest generation of wide field 
cameras on 6-10 meter class telescopes, near-IR obser- 
vations can now probe a large sample of young clusters 
with the sensitivity to detect objects well below the hy- 
drogen burning limit, the angular resolution to resolve 
high density groupings of stars, and the field of view 
necessary to observe the distribution of stars over multi- 
parsec distances. Furthermore, the Spitzer Space Tele- 
scope is now providing detailed images of young clusters 
in the mid-IR, allowing us for the first time to identify 
young stars with disks and infalling enve lopes efficiently 
in clusters out to 1 k pc an d beyond fe.g. . IMegeath et all 
2004; IWnitnev et alMOOl . 

Theoretically, the problem of star formation in clusters 
is perhaps an even greater challenge, as a realistic analy- 
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sis must include a network of non-linear processes, such 
as turbulent motions and shocks in the natal molecular 
cloud, potential collisions and multi-body interactions 
between young stars and protostars, and the destructive 
influence of radiation and winds from the newly formed 
stars. Numerical models of turbulent fragmentation have 
been successful at simulating the formation of clusters 
(or at least clusters of protostellar gravitationally bound 
cores) and these simulations are beginning to yield a 
wealth of theoretical predictions wh ich can now be in- 
vestig ated observationally fe.fi.. [Klessen fc Burkettl 

WM. IKlessen. Heitsch. fc Mac Low! 1200(1 

Bate Bonnell. fc Bromml 12003 iBonne 1. Bate, fc Vine) 
2003HGammie et alJ^OOa iLi et alJl2004D. For example 
numer ical simulations by iBate. Bonnell. fc Bromml 
( 2003) and IBonnell. Bate, fc Vine! l)2003j) of turbulent, 
initially spherically symmetric, Jeans unstable molecular 
clouds produce collapsing, non-uniform filamentary den- 
sity structures which are formed through the dissipation 
of kinetic energy by hydrodynamic shocks, allowing 
gravity to overcome turbulent support and force col- 
lapse. In these simulations, the gas eventually fragments 
into small (< 200 AU), self-gravitating, sub-virial 
clumps, which are interpreted as protostars. The models 
produce multiple groups of protostars which eventually 
merge to form large clusters. However, these simulations 
have obvious limitations: there are numerical problems 
associated with following the collapse of gas from parsec 
to 100 AU size scales, the models assume idealized initial 
conditions, and the models also do not take into account 
the feedback from the newly formed stars. In particular, 
the later stages of the evolution in these models may 
have questionable relevance to real clusters given that 
most young stellar clusters do not evolve into bound 
clusters, in part due to the disp ersal of most of th e 
binding mass before it forms stars ijLada fc L ada 2003). 

With increasingly sophisticated numerical models be- 
ginning to predict the spatial distribution and densities 
of young stars, it is an appropriate time to explore the 
configuration of young stellar clusters through wide-field 
infrared and submillimeter imaging, both to test current 
models and to guide the development of future mod- 
els. Through the analysis of the stellar surface den- 
sity in young (< 3 Myr) stellar clusters still associated 
with their natal molecular gas, and through a compari- 
son of the corresponding dust and gas density structure, 
we can study the initial configuration of young clusters 
and their subsequent dynamical evolution. We are now 
engaged in a multi-wavelength study of the structure of 
young stellar clusters, including ground-based near-IR 
data, Spitzer mid-IR imaging, and gro und-based submil- 
limeter and millimeter-wav e imaging ijGutermuth et ahl 
120041 iMeeeath et"aT]l2004|) . While testing cluster for- 
mation models is a long-term goal of the project, this 
paper is primarily focused on developing quantitative 
techniques for analyzing the spatial structure of young 
stellar clusters and applying these techniques to new ob- 
servations of three clusters. Subsequent papers will apply 
these techniques to the full sample of clusters observed in 
our survey. It is our hope that this work will encourage 
modelers to provide more accessible observable quanti- 
ties from their simulations. 

In this paper, we present an initial analysis of near-IR 
and submillimeter imaging of three clusters using near- 



IR images from FLAMINGOS on the 6.5 meter MMT 
telescope and 850 /im continuum maps from SCUBA on 
the JCMT. We restrict our near-IR analysis to an H — K 
color-based extinction analysis and the distribution of 
sources detected in the if -band. The if -band is th e 
most sensitive to embedded young stars i|Megeatbl il994V 
and thus gives the most unbiased view of the distribution 
of point sources. Although the membership of individ- 
ual sources cannot be established by the if-band data 
alone, statistical methods can be used to infer the overall 
distribution of young sources. This approach is particu- 
larly appropriate for studying the densest parts of young 
clusters, where the density of young stars dominates over 
the density of field stars in the line of sight. In a future 
paper, we will use data from Spitzer's Infrared Array 
Camera to find the distribution of sources with infrared 
excess emission in all three of the regions presented; this 
approach is particularly effective in probing regions with 
a lower density of young stars than the background stel- 
lar density. 

The three cluster regions chosen for this analysis ex- 
hibit signs of ongoing star-formation such as active 
molecular outflows, H2O masers, and associated IRAS 
point sources with significant far-IR luminosities ranging 
from 227 to 5.68 x 10 3 L Q (see TableUJ. A brief summary 
of the current knowledge of these three regions follows. 

1.1. GGD 12-15 

GGD 12-15 is a young cluster region associated 
with the Monoce ros molecular cloud, lo cated at a dis- 
tance of 830 pc ijHerbst fc Racing fl 9 76). It was origi- 
nally identified from several patches of optical reflection 
nebulosity with associated Herbig-Haro (HH ) objects 
l)Gvulbudaghian. Glushkov. fc Denisv ukl l 1978ft . There is 
a wide-spread (1.13 pc) stellar cluster detectable at near- 
IR wavelengths with a n estimated 134 objects (jHodaod 
1 19941 iCarpenteil |2000|) and a moderate velocity bipo- 
lar molecular outflow oriented northwest-southeast cen- 
tered on an H?Q maser and secondary peak in 800 /jm 
emission (iLittle. Heaton. fc DentHl990f) . A compact H II 
region, with flux consistent with excitation by a B0. 5 
zero-age main sequence star, was detected in radio con- 
tinuum emission l|Rodriguez et al.|ll980fl. The primary- 
peaks of 800 /im dust emission ({Little. Heaton. fc Dent 
1990), 13 CO emission, and C 18 Q emission |RkIge"ljtHd1 
2003) are also approximately centered on the H II re- 
gion. Furthermore, a cluster of point-like radio contin- 
uum sources has been detected in the vicinity of the H II 
region l)G6mez. Rodriguez, fc Garavll2000[ l20f)l . Some 
of these objects have near-IR and mid-IR spectral en- 
ergy distributions co nsistent with embedded protostars 
l|Persi fc Tapiall2003]) . 

1.2. IRAS 20050+2720 

IRAS 20050+2720 is a relatively small but very active 
site of embedded star formation originally identified as 
an IRAS point source. It is associ ated with the Cygnu s 
Rift, located at a distance of 700 pc (fWilking et al.11 989'1. 
Multipolar outflows have been found centered on the 
main IRAS source p osition, one of particularly high ve - 
locity (+/- 7 0km /s) ijBachiller. Fuente. fc Tafallal 1995ft . 
iChen et alJ (^97) performed a near-IR investigation of 
this object (completeness limit of K = 16.2), and re- 
ported an estimated 100 cluster objects arranged in three 
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tightly packed subclusters which they labeled A, B, and 
C. The most reddened of these is subcluster A, located 
at the IRAS point source p osition, as well as the peak of 
13 CO and C 18 emission l|Ridge et al.ll2003j) . The area 
around this subcluster has significant near-IR nebulos- 
ity, and two embedded sources were detected in L-band 
(3.5 /im) and narrow-band M (4.8 /im) imagi ng that were 
not d e tected at shorter n ear-IR wavelengths ijChen et al.l 
1997). iChini et alJ l)2001[) performed a submillimcter and 
millimeter study and found a strong emission source at 
the IRAS point source position, with a ridge extend- 
ing to the south. They also detected a diffuse ridge to 
the northwest at 1.3 mm and three point-like millimeter 
sources outside the 2.3' field of view of their 450 /im and 
850 /im maps. These are likely protostars or pre-stellar 
cores associated with the cluster. 

1.3. NGC 7129 

NGC 7129 is a large region o f reflection ne bulosity that 
is located 1 kpc from the Sun ijRacinell 1968(1 . It is illumi- 
nated by t wo Be stars and an associated cluster of lower 
mass stars ( HodaDD ll994l) . This cluster is located within 
a cavity west of a kidney-shape d molecular cloud w hich 
is sharply defined in b oth CO iRidee et alJ I2003T) an d 
submillimeter emission l)Font. Mitchell, fc SandelM 2001'). 
A nebulous filament which traces the cavity wall is the 
dominant feature of the region at both optical and near- 
IR wavelengths. The activity of forming young stel- 
lar objects has been detected in several places within 
the molecular cloud. Within the cavity wall is a third 
massive star, LkHa 234, and its embedde d companion 
l|Weintraub et alJ 119961 iCabrit et aD 11997ft A je t dis- 
covered at visible wavelengths llRav et alJ 11990ft w ith 
near-IR and mid-IR counterparts l|Cabrit et al]ll997ft is 
centered on these two objects. The jet is blue-shifted 
and points southwest into the cavity, suggesting that 
it may be the counterjet to the red-shifted molecu- 
lar outflow to the northeast. Another deeply embed- 
ded intermediate mass protostellar object, FIRS2, is lo- 
cat ed three arcminutes to the sou th of the main clus- 
ter l|Eiroa. Palacios. & Casalll99"8|) . FIRS2 is located at 
the re gion's primary peak in 13 CO emission (|Bechis et alJ 
11978ft and is associated with a multipolar molecular out- 
flow, sug gesting that it may have one or more you ng com- 
panions l|Fuente et al.l200HlMiskolczi et all2001ft . Many 
other Herbig-Haro objects have been detected around the 
regi on which are associated with sites of outflow activ- 
ity l|Edwards fc Sneljl983tlHartigan fc Ladall985ft . Sev- 
eral large-scale outfl ows have been observed in molecular 
hydrogen emission l|Eisloffelll200rl . Spitzer 3 to 8 /im 
and 24 /im imaging has revealed a significant extended 
cluster of embedded Class I and II young stellar ob- 
jects to the northe ast of LkHa 234 ijMuzerolle et al.l206"i 
iGutermuth et alJl2004ft . 

2. OBSERVATIONS 

Near-IR observations in the J (1.25 /im), H (1.65 /im), 
and K s (2.15 /mi) wavebands of NGC 7129 5 and 
IRAS 20050+2720 were obtained by the authors on 
June 15-16, 2001 using the FLAMINGOS 6 instrument 

5 The J and H band obs ervations of NGC 7129 w ere originally 
presented by the authors in Gutcrmuth ct al. (2004) 

6 We wish to acknowledge the superb instrumentation developed 
by Richard Elston, who helped transform the study of star forma- 



ijElstonl [1998) on the 6.5 meter MMT Telescope 7 . 
FLAMINGOS employs a 2048 x 2048 pixel Hawaii II 
HgCdTc detector array. The platescale of FLAMINGOS 
on the MMT is 0".167 pixel -1 , and thus the field of view 
is 5'. 7 x 5'. 7 per image. An overlapping 2x2 + 1 center 
position raster was used to make a 9'. 5 x 9'. 5 mosaic. 
Five dithered mosaics were obtained in each band with 
15 seconds exposure time per image, for a total exposure 
time of 75 seconds per band at J, H, and K s . Similar 
observations in the J, H, and K (2.18 /im) (instead of 
K s ) wavebands of GGD 12-15 were obtained on Febru- 
ary 23, 2002, also with FLAMINGOS on the MMT Tele- 
scope. The same raster pattern was used, but the expo- 
sure time per image was 30 seconds. We obtained eight 
dithered mosaics in J-band for a total of 240 seconds of 
exposure time, and four each in H and A'-band for a 
total of 120 seconds of exposure time in each band. Con- 
ditions were photometric for all observations, and stellar 
point spread function sizes ranged from 1" to 1.1" full 
width at half maximum (FWHM) in all bands. Only the 
H and K (Ks) data are presented here. 

Basic reduction of the near-IR data was performed 
using custom IDL routines, developed by the authors, 
which include modules for linearization, flat-field cre- 
ation and application, background frame creation and 
subtraction, distortion measurement and correction, and 
mosaicking. Point source detection and synthetic aper- 
ture photometry of all point sources were carried out 
using Phot Vis version 1.09, an IDL GUI-based photom- 
etry visualization tool developed by R. A. Gutermuth. 
Phot Vis utilizes DAOPHOT modules ported to IDL as 
part of the IDL Astronomy User's Library l|Landsmanl 
1993). By visual inspection, those detections that were 
identified as structured nebulosity were considered non- 
stellar and rejected. Radii of the apertures and inner 
and outer limits of the sky annuli were 1", 2", and 3.2" 
respectively. FLAMINGOS photometry was calibrated 
by minimizing residuals to corresponding 2MASS detec- 
tions, using only those objects with H — K s < 0.6 mag to 
minimize color differences in the 2MASS and FLAMIN- 
GOS filter sets. RMS scatter of the residuals between 
the two datasets ranged from 0.11 to 0.14 magnitude 
for the data presented. Median photometric uncertainty 
intrinsic to our FLAMINGOS measurements was 0.03 
magnitude, and our adopted maximum uncertainty tol- 
erance was 0.1 magnitude. Photometric magnitudes for 
stars which were brighter than the range of magnitudes 
corrected by the linearization module were replaced with 
2MASS photometry. The 90% completeness limits are 
derived by adding successively dimmer sets of synthetic 
stars created by using a Gaussian with a FWHM equal 
to that of the observed point sources in each mosaic, 
extracting their fluxes using the procedure described 
above, and rejecting objects with anomalous measure- 
ments. Measurements were considered anomalous if their 
photometric uncertainties, cr, were greater than 0.1 mag- 
nitude, or if the recovered magnitude was more than 3<7 
brighter than the input magnitude, indicating that the 
synthetic star was coincident with a brighter star in the 

tion during his too short life. 

7 Observations reported here were obtained at the MMT Ob- 
servatory, a joint facility of the Smithsonian Institution and the 
University of Arizona. 
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image. The magnitude at which 90% of the stars are 
recovered is recorded in Table 

The three clusters were observed with the SCUBA sub- 
millimeter camera at 850 on the James Clerk Maxwell 
Telescope 8 (JCMT) on Mauna Kea, Hawaii during a se- 
ries of runs from June to November 2003 in moderate 
weather conditions (r(850 /im) = 0.2 — 0.5). The data 
were obtained in scan-mapping mode whereby a map 
is built up by sweeping the telescope across the source 
while chopping the secondary at a series of throws, 30", 
44", and 68". To minimize striping artifacts, each set 
of three sweeps were performed alternately in right as- 
cension and declination. The SURF reduction package 
ijJenness. Lightfoot. fc Holland! IT998) was used to flat- 
field the data, remove bad pixels, and make images. Al- 
though effective at removing the high sky and instrument 
background, this observing technique necessarily loses all 
features on scales larger than the largest chopper throw, 
68". The calibration was performed by using the same 
observing technique to observe Uranus (when available) 
or standard sources CRL 618, CRL 2688, or IRAS 16293. 
These data are presented here for morphological compar- 
ison to stellar spatial distributions. A full analysis of the 
SCUBA maps will be presented in a future paper. 

3. METHODOLOGIES FOR STELLAR SURFACE DENSITY 
ANALYSIS 

Discriminating embedded young stars in clusters from 
field stars is difficult since the colors of the young stars 
are similar to those of field dwarfs and giants, and since 
reddening due to extinction dominates the colors of these 
sources. Young stars with strong IR-excesses from disks 
and envelopes can be identified, and we will adopt this 
approach by combining Spitzer and ground-based near- 
IR data in a forthcoming publication. However, this tech- 
nique cannot identify diskless young stars. An alterna- 
tive is to use statistical methods based on stellar surface 
densities to isolate regions on the sky that are dominated 
by cluster stars and measure several properties of the 
cluster stellar distribution. In the following subsections, 
we describe methods for isolating cluster-dominated re- 
gions of sky, modeling field star contamination there, and 
measuring several properties of the underlying stellar dis- 
tributions of the clusters. We apply these met hods to 
our observations of the three clusters. In Section l3~Tl we 
map the extinction through the molecular clouds associ- 
ated with each cluster, a key step in modeling the field 
star contamination. In Section [3.21 we use if-band mag- 
nitude histograms and the extinction maps to estimate 
the field star contamination and the number of cluster 
members in the high density central regions of each clus- 
ter. Azimuthally averaged stellar surface density profile 
fitting is used to characterize the cluster sizes and loca- 
tions in Section 13.31 The azimuthally averaged profiles 
prove to be incomplete characterizations of the spatial 
distributions. As shown by an analysis of the azimuthal 
distributions of stars in Section l3.4l there are significant 
deviations from circular symmetry in two of the three 
clusters. Finally, in Section [3.51 we present a method for 

8 The James Clerk Maxwell Telescope is operated by The Joint 
Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands Organ- 
isation for Scientific Research, and the National Research Council 
of Canada. 



mapping stellar surface density distributions using near- 
est neighbor distances, yielding another visualization of 
the often asymmetric structure present in the stellar dis- 
tributions of these clusters. 

3.1. Extinction Mapping 

Foreground and background stars are contaminants to 
studies of star clusters, obscuring the stellar distribu- 
tions of interest. Young cluster environments are rich 
with dust, reducing the density of background stellar 
contamination. However, the dust in these environments 
is often distributed non-uniformly, making it far more 
challenging to accurately characterize the distribution of 
detectable background stars. 

In order to correct for field star contamination, we 
must first map the distribution of extinction in our clus- 
ter regions. The SCUBA 850 /xm maps presented later 
in this paper cannot be used alone for this, as dust emis- 
sion which varies on scales larger than the largest chop- 
per throw (68") is lost due to the chopping sky sub- 
tracti on (see Section ?? ). Available CO emission line 
maps l|Ridge et al.ll2003^ are also not ideal for this pur- 
pose for many reasons, including low signal to noise in 
diffuse regions (C 18 0), optical depth issues leading to 
saturation in dense cores ( 13 CO), low spatial resolution 
(50" beam size in both maps), CO depletion due to 
chemistry or freeze-out in dense cores (both), or vary- 
ing excitation temperature (both). Instead, we derive 
the Ak maps for this work using the H — K colors of 
all stars detected in bot h bandpasses ijLada et alJl!994l: 
lLombardi fc Arvesll2001lK Similar to the stellar surface 
density mapping method described later on, we choose 
to consider the H — K values of the nearest N = 20 
stars at each point on a uniform 2" grid to match the 
SCUBA map sampling. Using an interative outlier re- 
jection algorithm, we compute the mean H — K value 
from the N nearest stars to that grid point. The algo- 
rithm calculates the mean and standard deviation of the 
H — K values, rejects any stars with an H — K value 
3cr from the mean, and then reiterates these steps until 
the mean converges. Outlier rejection in this application 
should primarily reject effectively unextinguished fore- 
ground stars, because the heliocentric distances to the 
clusters studied are fairly small. The resulting extinction 
map is convolved with a 15" Gaussian kernel to match 
the beam size of the SCUBA data and then converted 
to Ak using the reddening law of iCohen et al.l l)1981fl 
{Ak = 1.385((iJ — K) bs — (H — K) int )) and an assumed 
average intrinsic color (H — K)i nt — 0.2. 

This technique can underestimate the extinction since 
the average (H — K) color can include a contribution 
from young stars embedded in the molecular cloud; these 
stars have extinctions below the total line of sight ex- 
tinction through the cloud and bias the extinction to- 
ward lower values. This bias is particularly important 
in regions of high extinction, where the number of de- 
tected background stars is small. (Note that excess K- 
band emission from YSOs will not have a significant 
effect compared to scatter about the assumed intrin- 
sic H — K color.) To address this bias, the SCUBA 
maps are used to compare against the Ak map and 
to improve it where necessary. The SCUBA fluxes 
are multiplied by Ak / F^q — 0.272 mag pixel mJy -1 
(2" x 2" pixels) to convert to Ak and resampled onto the 
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same grid as the iJA-derived Ak map. The Ak/F S50 
conversion factor is derived assu ming Nn^/Ay = 
10 21 cm" 2 mag" 1 (see A ppendix D oflBertoldi fc McKed 
11992(1 . A K /A V = 0.090 ilCohen et al.lll981fl. grain mas s 
opacity k 8 so = 0.01 cm 2 g" 1 ((Pollack et alJ Il994fl . 
and dus t temperatures of T = 30 AT for all three 
clusters (iLittle. Heaton. fc Dendll99Ti iChini et al.ll200H 
IFont. Mitchell, fc Sandel]ll2001lh If any SCUBA-derived 
Ax measurements are more than 30% above the corre- 
sponding _ff A-derived Ak values, the higher SCUBA- 
derived value is used. This hybrid Ak map is used for 
the background contamination simulation. The original 
i? A-derived Ak map is preserved for use in a simula- 
tion of the clusters themselves, described in Section EQ1 
below. Both maps are expanded and resampled using bi- 
linear interpolation to match the near-IR data in spatial 
scale and position on the sky. See Figs. ^ d &l21f or the 
final extinction maps for GGD 12-15, IRAS 20050+2720, 
and NGC 7129, respectively. 

3.2. K-band Apparent Magnitude Histograms 

A-band Apparent Magnitude Histograms (KMH) 9 are 
commonly used tools in the study of young stellar clus- 
ters. In the literature, these histograms are commonly 
referred to as A"-band luminosity functions; for clus- 
ters with a known distance and negligible intracluster 
extinction, each apparent magnitude bin corresponds to 
a unique luminosity range in the observed wavelength 
band. The A-band is typically used since it is the least 
affected by extinction and the most sensitive to reddened 
young stellar objects. Thus, the KMH has been used to 
study the initial mass function of young stellar clusters. 
In this paper, we limit our analysis first to determin- 
ing the magnitude at which the contamination from field 
stars begins to dominate the A-band number counts, and 
second, to estimating the total number of stars in each 
cluster. 

We first define two regions in the observed field, the 
cluster-dominated region of the field and an off-field, dis- 
placed from the cluster, where we can obtain an estimate 
of the amount of field star contamination. We assume 
that the stars in the off-field are representative of the 
field stars over the entire observed field. The cluster- 
dominated region is defined via an analysis of the radial 
density profile (see Section 13.3(1 . Since the determina- 
tion of the radial density profile also depends on a cutoff 
magnitude at which the field stars begin to dominate the 
number counts, several iterations were needed to con- 
verge to a consistent magnitude cutoff and cluster region. 
The off-field is selected based on relative lack of CO emis- 
sion ((Ridge et al.l l2003). lack of submillimeter emission, 
lack of significant or structured extinction in the hybrid 
Ak maps derived in Section 1X11 ab ove . and overall lack of 
stellar surface density structure (see Section l3~K(l . Sepa- 
rate KMHs of the stars detected in both regions are then 
constructed using a bin size of 0.5 magnitudes. 

To determine the magnitude at which the field star 
contamination begins to dominate the number counts, 
we use the observed KMH of the off-field to model the 
magnitude distribution of the stellar contamination in 

9 To simplify the manuscript, all references to the FLAMINGOS 
if -band and if s -band Apparent Magnitude Histograms will be 
referred to as KMH regardless of the specific filter used. 



the cluster region. We assume that there are no cluster 
members in the off-field. To account for the higher ex- 
tinction toward the cluster, the background stars must 
be artificially extinguished. We separate the foreground 
and background components of the off-field KMH by 
estimating the KMH o f the foreground stars using the 
iWainscoat et alJ 1(1993(1 galactic star counts model and 
subtracting it from the observed off-field KMH. We then 
measure the difference between the mean cluster region 
and off- field extinction, SAk, using the hybrid Ak maps 
(Section 13.1(1 . This difference is applied by shifting the 
background KMH by SAk- Then the foreground contam- 
ination model is added to the shifted background model, 
and the total is scaled by the ratio of the areas of the 
cluster region and the off-field. The observation-derived 
stellar contamination estimates for the clusters presented 
agree to within statistical errors with t hose derived from 
the g alactic contamination model of IWainscoat et alJ 
(1993), once its background component has been shifted 
by the average Ak of the cluster region. 

The resulting stellar contamination model histogram 
is subtracted from the raw cluster KMH. By inspection 
of the resulting difference histogram, we note the magni- 
tude at which the cluster region becomes dominated by 
field star contamination. Specifically, we use the faintest 
magnitude bin to have greater than 3o~bkg stars after field 
star subtraction, where o~bkg is derived by assuming that 
Poisson counting statistics holds for the number of stars 
in each bin of the raw off-field KMH. We adopt this mag- 
nitude as the limiting magnitude for our data if it is 
brighter than the 90% completeness limit, otherwise the 
end of the first complete magnitude bin brighter than the 
90% completeness limit is used. Stars fainter than the 
limiting magnitude are not included in the stellar surface 
density analyses presented in this work. 

For the purposes of determining the total number of 
stars and studying the spatial distribution of stars in 
each cluster, a more sophisticated treatment of the field 
star contamination is required which takes into account 
the spatially varying extinction. To estimate the po- 
sition dependent magnitude distribution of field stars, 
we construct Monte Carlo simulations of the foreground 
and background stellar distributions in each of the re- 
gions observed. We use the KMH measured from the 
off-field, minus the expected for eground KMH from the 
model of Wainscoa t et al.l ((1993) , to determine the num- 
ber of background stars per magnitude per square pixel 
of the image plane. These are assumed to have mean 
extinction Ak,o, measured directly from the off-field of 
the hybrid Ak map, as before. We then uniformly popu- 
late an image plane equal in angular size to our observa- 
tions with a Poisson-deviate number of stars with mean 
equal to the total of the off-field background KMH, ap- 
propriately scaled to the area of the image plane. The 
magnitude of each star is randomly sampled from the 
off-field KMH and the appropriate positionally depen- 
dent extinction, SAk — Ak[x,v] — Ak,0i is applied to 
each star. All stars dimmer than the adopted magni- 
tude limit for the given cluster observations are rejected. 
To this sample we add a similar uniformly distributed 
population of unexting uished foreground stars, using the 
IWainscoat et alJ {Ti) 9."> ■ model foreground KMH to deter- 
mine their density and magnitudes. We perform Nqm 
iterations of the contamination model, where we have 
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chosen to use Ncm — 1000 for this work. These itera- 
tions are combined for use in the final field star contam- 
ination characterizations which will be used throughout 
the rest of this work. Depending on the dust distribu- 
tion, these simulations can predict a higher amount of 
contamination than that obtained using an average Ak 
(for example, see the middle plot of Fig. 0J| 

The final contamination histogram is measured di- 
rectly from the stars that fall in the cluster region of 
our simulated field star distributions, and divided by 
Ncm- This is subtracted from the raw cluster KMH 
we observed. The star counts in Table EH are calculated 
by totaling the contamination-subtracted cluster region 
histogram for all magnitudes brighter than the limiting 
magnitude adopted (see the final plot in Fig.0J). 

3.3. Radial Density Profiles 

Many studies of young stellar clusters in the litera- 
ture have analyzed azimutha lly averaged radial pro files of 
stellar surface density (e.g.. lCarnenter et al.l ll997'). This 
method allows the resulting trend to be fit by a wide 
variety of functions and naturally yields a size estimate 
for the cluster. Unfortunately, physical interpretations 
derived from these analyses are often suspect because 
many young clusters a ppear to exhibit non -uniform and 
asymmetric structure l)Lada fc Ladal 12003^ . As we will 
show, azimuthal averaging of elongated clusters can lead 
to deceptively smooth, well behaved radial density pro- 
files. Thus for this study, we limit our radial density 
profile analyses to estimation of the size and location of 
the cluster core and cluster-dominated region. 

To generate azimuthally averaged radial stellar surface 
density profiles for each cluster, we first choose a center 
point that falls approximately at the center of the pro- 
jected extent of the stellar cluster. This choice is then 
refined by taking the mean of the Right Ascension and 
Declination of all stars detected within one arcminute of 
the initially chosen center point. Based on this refined 
center point, we count the number of stars in concentric 
annuli of equal radial extent and divide by the area of 
each annulus. The radial bin size used in this work is 
equivalent to 0.034 pc at the assumed cluster distances. 
It was chosen to achieve adequate spatial resolution while 
maintaining acceptable statistical weight for the inner 
stellar density measurements. Given N stars in the i th 
annulus, the stellar surface density for that radial bin is 
computed as follows: 

Ni 



We fit the measured radial surface density histogram, 
weighting each bin by the number of stars it contains, 
using a function of the form: 

cr(r) =a + a ie { - r/ra) 

This function implies that field star contamination has 
a uniform surface density oq. As mentioned in Sec- 
tion 13.11 background stellar densities will vary signifi- 
cantly in regions of variable extinction. Using the stel- 
lar distributions from the Monte Carlo contamination 
models produced in Section [3.21 we measure azimuthally 
averaged radial profiles from the combined Ncm stel- 
lar contamination models using the identical bin sizes 



and center points used for the observations. The re- 
sulting densities are then divided by Ncm, and sub- 
tracted from the corresponding observation-derived den- 
sities. These contamination-subtracted profiles are then 
fit with the above function, ensuring that the cluster pro- 
files are properly characterized independently of nonuni- 
form background contamination (cf. Fig. [SJ . 

We define the cluster-dominated region as a circle cen- 
tered at the refined center point with radius of the point 
at which the contamination-subtracted profile fit de- 
scends below 3(7 above the average residual field star 
density (see Table EU, where a is computed directly from 
the outermost six radial bins. This isolates the region of 
elevated stellar density, and thus the extent of the region 
that is dominated by cluster members. We note empir- 
ically that three times the e- folding length, ro, approxi- 
mates this distance well in the two centrally condensed 
clusters presented, but not in the case of a more diffuse 
cluster. This value is also noted in Table 01 as the cluster 
core radius. 

3.4. Azimuthal Distribution Histograms 

Radial density profiles and azimuthally averaged dis- 
tributions of the stellar density in clusters do not take 
into account the asymmetric and often clu mpy structure 
obser ved in many young stellar clusters l)Lada fc Ladal 
2003). Indeed, such methods actually erase asymmetric 
structure which may yield vital insight into the process of 
clustered star formation. To examine the azimuthal dis- 
tribution of sources and determine whether the clusters 
presented show a high degree of azimuthal asymmetry, 
we construct azimuthal distribution histograms for each 
of the three clusters. Specifically, we count the number 
of stars in the cluster-dominated region that are located 
in 45° position angle bins sampled at 22.5°, the Nyquist 
rate. 

To quantify the azimuthal asymmetry of a cluster, we 
first compute the standard deviation of the histogram 
about the mean number of stars per azimuthal bin. We 
assume that this is the combined statistical noise of the 
cluster distribution and the field star distribution. We 
define the azimuthal asymmetry parameter (AAP) of a 
young stellar cluster as the ratio of the measured stan- 
dard deviation of the azimuthal distribution to the square 
root of the mean number of stars per azimuthal bin. 
Given m uniform overlapping (360/(m/2))° azimuthal 
bins with m stars in the i th bin and mean \n\ = N / (m/2) 
stars per bin, the raw AAP is computed as follows: 



m (rii — \n\) 2 
m — 1 



Em 
i=l 



AAPn 



A circularly symmetric cluster has a random azimuthal 
distribution. Therefore its azimuthal distribution his- 
togram should be uniform to within Poisson noise, as 
should a uniform distribution of contaminating field 
stars. This implies that a combination of both also has 
a uniform azimuthal distribution. Given successful iso- 
lation of a subset of stars in a region which is dominated 
by cluster members, either via the statistical methods 
described above or via true population restriction such 
as the consideration of stars exhibiting excess infrared 
emission only, a quantitative statistical argument for the 
presence of significant asymmetric structure can be made 
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using the measured AAP. Note that better isolation of 
cluster members reduces additional noise from the ran- 
dom field star distribution, yielding a stronger statistical 
argument. 

One additional complication is that variable extinc- 
tion can introduce asymmetry in the background stel- 
lar distribution. This must be accounted for in order 
to properly isolate any asymmetry that is due to the 
cluster distribution. Similarly to Section 13.31 we can 
measure the distribution of stellar contamination in az- 
imuth using the Monte Carlo simulations performed in 
Section l3~2l We measure the azimuthal distribution his- 
togram of the stars that fall in the cluster region bound- 
ary of the Ncm contamination models, and divide each 
bin by Ncm- Subtracting this mean model distribution 
from the observed one removes the systematic asymme- 
try expected in the background stellar distribution, but 
lacks the appropriate Poisson scatter. Thus it is impor- 
tant to note that noise from the field star distribution is 
not removed by this method. Given hi stars in the i th bin 
of the field star model azimuthal distribution and mean 
field stars per bin of \b\, the corrected AAP is the mea- 
sured standard deviation of this subtracted distribution 
over the square root of the mean of the unsubtracted dis- 
tribution, \n\, plus a small contribution from the Poissan 
noise of the mean field star model itself: 



AAP S „ 



EI 



{rii — bi — \n— b\ ) 2 
= 1 m-1 



y/\n\ + \b\/N CM 



In order to argue for the presence of significant asym- 
metric structure, we must detect azimuthal stellar sur- 
face density variations larger than that which is likely via 
Poisson counting statistics alone. To determine the prob- 
ability that a measured AAP is consistent with a ran- 
dom azimuthal distribution, we have performed Monte 
Carlo simulations of 10 6 young cluster fields for each of 
the regions presented. The foreground and background 
components of each simulated field are generated identi- 
cally to our field star contamination models, as described 
in Section 13.21 In addition, the clusters themselves are 
also modeled. For the reference magnitude distribution 
for the clusters, we have chosen to use the well-sampled 
.ftT-band luminosity function (KLF ) of th e Orion Nebula 
Cluster (ONC) from lMuench et aH p0021. The exponen- 
tial fits of the subtracted radial profiles (see Section l3~3)l 
are used as the reference radial density distributions for 
our model circularly symmetric clusters. The number of 
cluster members in a given model is a Poisson-deviate of 
the integrated exponential component of the radial den- 
sity profile fit. The stars should be distributed uniformly 
in position angle about the measured cluster center point 
of the observations, but we must account for reduced ef- 
fective sensitivity to low mass cluster members in higher 
extinction areas. To do this, we compute an azimuthal 
distribution weighting function for a given radius by inte- 
grating the assumed KLF up to the adopted magnitude 
limit minus the HK-derived Ax value at that radius 
over all position angles. The normalized results yield the 
position angle probability distribution functions for the 
range of radii possible in our field. 

Each of the 10 6 young cluster fields generated is ana- 
lyzed identically to the observations. The refined clus- 
ter center points are computed, using the center point 



measured from the observations as the initially chosen 
point which is then further refined. Radial profiles are 
measured and the contamination models used to treat 
the original observations are also used to compute and 
remove radial variations in stellar contamination. The 
resulting radial profile is fit as described in Section 13.31 
with those simulations with nonconverging fits rejected 
as unphysical and regenerated from the model 10 . The 
cluster-dominated region is thus defined as in Section ^. 31 
The azimuthal distribution is then measured, and we use 
the same mean of Ncm stellar contamination models, 
recentered on that field's refined cluster center point, to 
remove any systematic variation due to structured ex- 
tinction. The AAP of each iteration is then computed, 
compiled into a histogram, and then normalized to yield 
the probability distribution of the AAP over all the sim- 
ulations. 

From these distributions, we can determine the prob- 
ability that the measured AAP of the cluster of stars 
observed is consistent with a random azimuthal distri- 
bution, and is thus an approximately circularly symmet- 
ric cluster. A high measured AAP suggests a very small 
probability that the cluster is circularly symmetric, im- 
plying the presence of significant asymmetric structure. 
As an example, inspection of the GGD 12-15 proba- 
bility distribution in Fig. suggests that a measured 
AAP of 1.588 has a -0.008% probability of being con- 
sistent with a circularly symmetric distribution. Thus 
the presence of azimuthal asymmetry can be argued for 
clusters with measured AAP values significantly greater 
than |A<4P S j m |, while those with AAP values close to 
|j4AP s j m | are consistent with a circularly symmetric dis- 
tribution. The resulting field star model subtracted 
lAAPsiml values from the simulations for GGD 12-15, 
IRAS 20050+2720, and NGC 7129 were 0.859, 0.869, and 
0.894 respectively. 

3.5. Stellar Surface Density Maps 

Ultimately, a stellar density distribution mapping 
method must be employed to investigate the relation- 
ship between forming stars and their natal molecular 
gas distribution. Given the wide range of stellar sur- 
face densities in a young cluster region, a method that 
employs adaptive smoothing lengths is ideal for achiev- 
ing both high dynamic range density measurements and 
high spatial resolution in locations of high stellar surface 
density. For this analysis, we have ch osen to use a varia- 
tion o n the nearest neighbor method ijChristonher et all 
1998) to construct the stellar surface density maps for 
each cluster. This method exhibits many of the benefits 
of the GATHER algorithm outlined by iGladwin et alJ 
(1999), but lacks weighted averaging of multiple smooth- 
ing lengths. At each sample position in a uniform 
grid we measure r^(i,j), the projected radial distance to 

10 Nonconvergent profile fits most often occur in clusters with 
low ratios of peak surface density to field surface density, such 
as NGC 7129 in this work. Since the resultant radial profile of a 
given iteration of the model is allowed to vary according to Poisson 
scatter, clusters with less central condensation may produce stellar 
distributions that are indistiguishable from flat density plateaus 
or other non-exponential geometries. Model stellar distributions 
that take on such profiles which lack central peaks and exponential 
character are unlikely to be identifiable as cohesive clusters distin- 
guishable from field stars by our stellar surface density isolation 
methods, and therefore we reject those iterations. 
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the N nearest star. For convenience, we chose a grid 
spacing of 3" for all stellar surface density maps pre- 
sented in this work. The local stellar surface density at 
each grid position is computed as follows: 



In Figure we present the inner 5' of the stellar den- 
sity maps derived with this method, using N = 5 and 
N = 10 (left and center columns respectively) to demon- 
strate how the maps change as a function of N. The 
most obvious difference between the two is loss of spa- 
tial resolution with larger TV. A more subtle difference 
is the total loss of some small, high density features in 
the N — 10 maps. Any high density subgroups that 
contain less than N stars are often lost or significantly 
diminished with this method, thus it is critical to choose 
N to be small enough to allow small, high density sub- 
groups to be properly represented in the map. Further- 
more, N must be large enough that the measured density 
structure is not dominated by binary or triple systems or 
random coincidences. For these reasons, we have chosen 
to use maps derived for N = 5 for this study. Peak 
surface density entries in Table |31 were obtained directly 
from the N = 5 maps. For independent comparison, the 
right column of Fig.EJhas equivalent grid resolution stel- 
lar density maps derive d using Gaussian kernel smooth- 
ing l)Gomez et al.l lT993') with a smoothing length of 8". 
Note that the density distribution morphologies from the 
N = 5 maps are discernable using this commonly used 
alternative method. 

Approximate volume densities can be derived from sur- 
face density measurements only by making assumptions 
about the distribution of stars in the line of sight. Since 
our definition of a cluster-dominated region is spatially 
confined to a circle, one might naively assume that the 
stars were distributed evenly along the line of sight over 
a length of the order of the cluster core radius. How- 
ever, toward peaks in the observed surface density, it is 
unlikely that the stars are distributed uniformly over the 
cluster diameter along the line of sight, and this den- 
sity would be a lower limit to the peak volume density. 
An alternative method is to derive volume density (see 
TableOU using the nearest neighbor distance rpj(i,j) de- 
fined above, thus assuming local spherical symmetry at 
those sample positions. We use N = 5 for these mea- 
surements. For grid position the nearest neighbor 
volume density is computed as follows: 



This method would overestimate the volume density in 
clusters where the stars are uniformly distributed over 
the line of sight, but would provide a more accurate mea- 
surement of the volume density toward observed surface 
density peaks resulting from a sharply peaked or clumpy 
volume density field. The method is particularly appli- 
cable when the number of stars in a volume density peak 
is larger than the number of surrounding clusters stars in 
the line of sight. This method is also more applicable in 
highly elongated or filamentary clusters. To minimize the 
effect of field star contamination, we apply this method 
only in regions were the surface density of cluster mem- 
bers is much larger than the density of field stars. The 



reported densities are only densities of detected stars; 
stars too reddened or too dim to detect in our data are 
not taken into account. 

4. RESULTS 
4.1. GGD 12-15 

By inspection of the contamination-subtracted KMH 
of the cluster region of GGD 12-15 (Fig.0J it is clear that 
the cluster population dominates field star counts at all 
magnitudes throughout the sensitivity range of our data. 
Hence we choose to use a limiting magnitude equal to the 
end of the first complete magnitude bin brighter than the 
90% completeness limit, K — 18. The contamination- 
subtracted radial density profile fit (Fig. EJ) yields a clus- 
ter core radius of 0.24 pc, but the large variation in the 
azimuthal distribution histogram (Fig. [SJ suggests signif- 
icant asymmetric structure in the cluster core. Indeed, 
there are two peaks in the azimuthal distribution which 
are separated by 180°, a signature of linear structure. 
This is confirmed by inspection of the stellar distribution 
(see Figs. |H1 and EJ), the stellar density map fFig.HUI). and 
the stellar density map slice plots (Fig. Illfl . Finally, the 
contamination-subtracted AAP of GGD 12-15 is 1.588, 
3.77cr above the simulated mean AAP, suggesting only a 
0.008% probability that the radially averaged azimuthal 
distribution observed is consistent with a circularly sym- 
metric configuration. Clearly, an azimuthally averaged 
profile is a poor characterization for the highly struc- 
tured nature of this cluster. 

The refined cluster core center point approximately co- 
incides with the location of peak density of the stellar 
density distribution map of GGD 12-15, showing that 
there is a well-defined and cohesive stellar surface den- 
sity peak at the center of the cluster core, even though 
there is significant structure evident in the stellar density 
map. The position of peak stellar density occurs between 
two peaks in the core dust emission. This suggests cither 
that the central stars are in the process of opening a cav- 
ity in the cloud core or that local extinction has in fact 
drastically reduced our sensitivity to stars in the areas of 
high dust density. Furthermore, we report the discovery 
of 850 fim dust emission filaments (see Fig.|5J) extending 
from the main core which has been p reviously investi- 
gated by EES^^iilS^^^Si^ (El9(j). The orientation 
of these filaments is nearly identical to the orientation of 
the elongated stellar distribution (see Fig. HHjl . We argue 
that this is evidence that star formation in GGD 12-15 
is primarily occuring along the molecular cloud filament 
and that the stars have not had adequate time to migrate 
significantly from their birth sites. 

4.2. IRAS 20050+2720 

The contamination-subtracted KMH of the cluster re- 
gion of the IRAS 20050+2720 field (Fig. Q2J shows that 
the cluster population becomes dominated by field star 
contamination for sources fainter than K s = 16, hence 
we adopt that as our limiting magnitude for this work. 
The radial density profile fit (Fig. ffify yields a cluster 
core radius of 0.24 pc, but as with GGD 12-15, there is 
significant evidence to suggest that this is a very poor 
characterization due to significant deviation from cir- 
cular symmetry (see Figs. El an d El an d note that 
iChen et all 119971 reported clear evidence of three dis- 
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tinct subclusters in IRAS 20050+2720, which they des- 
ignated subclusters A, B, and C). The refined cluster 
core center point of IRAS 20050+2720 is not located at 
the position of peak density in the stellar density map, 
but is instead located approximately 20" to the south- 
west. The azimuthal distribution histogram (Fig, 114(1 has 
an obvious peak corresponding to the contribution from 
subcluster A and a less prominent one corresponding to 
subcluster C (see Figs. 1171 and 11811 . The contamination- 
subtracted AAP of IRAS 20050+2720 is 1.421, 2.82cr 
above the simulated mean AAP, suggesting a 0.2% prob- 
ability that the azimuthal distribution observed is consis- 
tent with a circularly symmetric confiuration. It is clear 
that extended stellar distributions with subclustering to 
this degree are poorly described using methods that uti- 
lize azimuthal averaging. 

While the locations of peak stellar surface density and 
the refined cluster core center point are not coincident, 
both positi o ns ar e within subcluster B as defined by 
iChen et alJ l|1997tK Subcluster B is in a region lacking 
850 /im emission within an otherwise apparently colin- 
car filamentary structure in the dust emission map (see 
Fig.EJ- This subcluster has significantly less reddening, 
dispersion in H — K, and associated if -band nebulosit y 
as compared to subclusters A and C l|Chen et alJll997j) . 
which are both located in areas of significant 850 /urn 
emission. However, the stellar surface density map (see 
Fig. I17fl does not show a clear boundary between sub- 
clusters A and B, leaving open the possibility that these 
two subclusters are not distinct structures, but instead 
may be subregions of a single subcluster that have very 
different amounts of associated dust in the line of sight. 

Our wide-field 850 /im map is nearly iden tical in mor- 
pholog y to the 1.3 mm map presented in iChini et"aT1 
(|2001IK including the small cores named MM2 and MM3 
located south of the cluster region and the small un- 
named core to the northwest of the cluster. Given the 
lack of dust emission in the vicinity of subcluster B and 
the multiple high velocity molecular outflows detected 
in this region, our analysis suggests that it has very re- 
cently dispersed its natal molecular gas while star forma- 
tion proceeds along the rest of t he filament, in agre ement 
with the conclusions reached bv lChen et alJ 1(19971) . 

4.3. NGC 7129 

By inspection of the contamination-subtracted KMH 
of the cluster region of NGC 7129 (Fig. |HJ) we ar- 
gue that the cluster population dominates field star 
counts at most magnitudes throughout the sensitivity 
range of our data. Hence we choose to use a limiting 
magnitude equal to the end of the first complete mag- 
nitude bin brighter than the 90% completeness limit, 
K = 17.5. The contamination-subtracted radial den- 
sity profile fit (Fig. I20f) yields a cluster core radius of 
0.65 pc, much larger than the that of the other two clus- 
ters presented in this work. Also, the azimuthal distri- 
bution histogram (Fig. 12 1|) is quite featureless by com- 
parison to the other two clusters presented, suggesting 
that the cluster is approximately circularly symmetric. 
This can also be seen by inspection of the stellar distri- 
bution fFig. 122(1 which is clearly located primarily in the 
cavity defined by the 850 /im emission (see Fig. 123(1 . Fi- 
nally, the contamination-subtracted AAP of NGC 7129 
is 1.207, 1.57cr above the simulated mean AAP, suggest- 



ing a 5.8% probability that the azimuthal distribution 
observed is consistent with a circularly symmetric con- 
figuration. While this result suggests some asymmetry 
may be present, inconsistency with a circularly symmet- 
ric distribution is not statistically significant. 

The refined cluster core center point is not located at 
the position of peak density in the stellar surface density 
map (see Fig. I24[l . but is instead located approximately 
30" to the west centered on another local density maxi- 
mum. There are three additional local density maxima 
within the cluster region boundary, and all five have sim- 
ilar peak stellar surface densities (see Fig. |2H| ■ Clearly, 
there is no dominant central high stellar density core in 
NGC 7129. 

5. DISCUSSION 
5.1. The Initial Configuration and Evolution of Clusters 

As we have shown, characterizing young stellar clus- 
ters with azimuthal averaging methods obscures a key 
component in their analysis. Those clusters pre- 
sented that are still highly embedded, GGD 12-15 and 
IRAS 20050+2720, have a high degree of asymmetry in 
their stellar distributions compared to NGC 7129. Given 
the similarity of the asymmetric stellar distributions to 
the filamentary structures detected in 850 /im emission, 
we argue that the stellar structure observed is a direct re- 
sult of asymmetry in the distribution of the natal gas. In 
contrast, the NGC 7129 cluster is approximately circu- 
larly symmetric, and occupies a cavity in its natal cloud. 
Low average and peak stellar densities within the clus- 
ter, overall lack of significant stellar density structure, 
and a significantly larger cluster core radius are circum- 
stantial evidence that recent expulsion of the bulk of the 
gas mass from the cavity has allowed the cluster to dy- 
namically expand to its current state. 

In the two embedded clusters there is also evidence 
that the star-forming gas and dust are being disrupted. 
While stellar density morphology seems related to dust 
emission morphology in the embedded clusters, peak stel- 
lar densities are often anti-correlated to dust emission 
maxima, suggesting that either local extinction variation 
severely affects our sensitivity to stars or stellar feedback 
mechanisms are already dispersing natal material locally 
near high stellar densities. For example, in GGD 12-15, 
the peak of stellar density is located at a local minimum 
in 850 fj.rn emission, suggesting that the central stars may 
be opening a cavity in the dense central cloud core. Sim- 
ilarly, subcluster B of IRAS 20050+2720 is just outside 
the region of peak 850 /im emission, possibly indicating 
th at it too has dispe rsed much of its gas, as suggested 
bv IChen et alJ <|1997l) . Subcluster B may be entering a 
phase of dynamical expansion similar to what we sug- 
gest to explain the current state of the NGC 7129 clus- 
ter, while the remaining filament is still actively forming 
stars in subclusters A and C a nd the small mill imeter 
emission cores detected nearby l|Chini et al.1120 01'). 

Analytical investigations and numerical simulations 
of the dynamics of young clusters following rapid ex- 
pulsion of their natal gas have shown that dynam- 
ical expansion and dissolution can occur on short 
time scales in young clusters l| A dams fc Mversl 1200 It 
iKroupa. Petr. fc McCaughreanlll999l) . Indeed, over the 
course of 1 Myr an unbound star moving at 1 km/s can 
travel 1 pc, suggesting that the high stellar densities ob- 
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served in young clusters can only last while the system is 
gravitationally bound by the mass of the natal molecular 
cloud. The observations presented in this work further 
these arguments. Understanding the impact of forming 
stars on their parent cloud is crucial to understanding 
the complex dynamics of embedded cluster evolution. By 
analyzing stellar density distribution morphology in rela- 
tion to molecular cloud structure, observational analyses 
can more adequately address the clear link between star 
formation, gas expulsion, and the dynamics of the clus- 
ters and how these processes guide the evolution of young 
clusters. 

5.2. The Impact of Clustering on Disk and Envelope 
Evolution 

Simulations of massive protostellar disk interactions 
suggest that truncation and fragmentation of the disks 
are likely for close approach es with other stars on the or- 
der of the disk radius fe.g.. lBoffin et al1 ll998). A recent 
investigation into the effects of close approaches on classi- 
cal T Tauri disks as a function of a variety of parameters 
suggests tidal truncation an d mass loss are poss ible, as 
well as enhanced accretion (Pf alzner et al J 12005). Fur- 
thermore, the gravitational instability planet formation 
scenario may require a close approach with another star 
to perturb the protoplanetary disk to begin the process 
ijBossI 12002]) . If close approaches are likely for a signif- 
icant fraction of pre-main sequence stars in an embed- 
ded cluster, this would have a profound effect on disk 
evolution and perhaps planet formation in a clustered 
star-forming environment. 

Observations other than those presented here provide 
further clues that envelope and disk interactions may 
happen frequently and indeed may be significant in the 
evolution of disks and the formation of planets. Several 
disks around pre-main sequence stars in the Orion Neb- 
ula Cluster (ONC) have been di rectly observed with HST 
ijMcCauehrean fc O'Delll I1996|) : these disks have radii 
ranging from as small as 50 AU to as large as 1000 AU, 
and their central stars have derived ages of 0.8-3.0 Myr. 
Given the high stellar densities measured in the ONC, 
as well as the wide range of disk radii and sharp disk 
edges observed, tidal tru ncations from close encounter s 
are a likely explanation l)McCaughrean fc Q'Delll fl996'l. 
A 10 AU gap in the disk of ~1 Myr old T Tauri star CoKu 
Tau/4 has been discovered in observations from the 
Spitzer Infra red Spectrograph, suggesting the presenc e 
of a planet l|Forrest et alJ l200l iD'Alessio et"alll200l . 
This suggests that giant planet formation may occur very 
early in a protoplanetary disk's lifetime, an argument 
in favor of planet formation via gravitational instability. 
Current orbital models for the recently discovered dis- 
tant minor planet Sedna |Brown. Truiillo. fc Rabinowitd 
l2004j) argue for a close interaction from a passing star, 
suggesting that our solar system may have formed in a 
clustered environment. 

One of the strengths of the stellar surface density map- 
ping method adopted for this work is its ability to probe 
stellar densities over a wide dynamic range with ade- 
quate spatial resolution. While stellar collisions are very 
unlikely at the stellar densities reported in this work (cf. 
iBonnell fc Batell2002t and note that lower average stel- 
lar masses further increase the timescales in Figure 1 of 
that study) , close approaches on the order of protostel- 



lar envelope or large T Tauri disk radii (R = 10 3 A U, 
IMotte fc Andrei EFioTl IMcCa,ughrean fc 0'DelillT99fil) at 

these densities may be much more common. The time- 
scale for a given star to experience a close encounter when 
among a cluster of stars of stellar density n and veloc- 
ity dispersion Vdisp can be roughly approximated via a 
simple analysis outlined in iBinnev fc Tremaind <|1987|> . 
They express this time-scale, t co u, as: 

1 r- GM X 



- = I&^/ttuv disp Rl{\ 

tr.nll 



ID > 



All stars are assumed to have the same mass, M*, and the 
same interaction radius, i?*. For this analysis, we assume 
ill* = 0.5 Mq as an appropriate median stellar mass from 
young cluster initial mass functions reported in the lit - 
erature (e.g., [Muench et al. 2002; Luhm an et alJ 2003) . 
and Vdisp = 1 km s _1 based on Maxwellian velocity dis- 
persions derived from C 18 linewidths for these clusters 
(see Table nj. Figure l2"o1 shows the derived collisional 
timescales as a function of stellar densit y for two choices 
of R* , 10 3 AU for protostellar envelopes ijMotte fc Andrei 
I2001() and large T Tauri disks and 10 2 AU for classi- 
cal T Tauri disks (jMcCaughrean fc O'Delll 119961. The 
timescales considered must be less than or on the order 
of 10 5 to 10 6 yr to allow for a significant frequency of 
interactions to occur over t he length of a forming star's 
protostellar collapse phase ijKenvon fc HartmannllT995]) 
or the approximate length of a young c luster's active 
star-forming phase ijPalla fc Stahleill2000|) . respectively. 
For example, at a density of 10 4 pc -3 , protostars and 
T Tauri stars with large disks are likely to have a close 
encounter on the order of their envelope or disk radius 
on a timescale of ~10 5 yr, while classical T Tauri disk 
radius interactions occur much less often, on a timescale 
of ^10 7 yr. The latter timescale suggests that ~10% 
of stars at densities of 10 4 pc~ 3 would have a close en- 
counter on the order of 10 2 AU in ~1 Myr, if the high 
stellar density environment can indeed last that long. 

Utilizing our stellar density maps and cluster member- 
ship and field star contamination analyses, we can esti- 
mate the prevailing stellar density environment around 
the stars in the cluster regions relative to the interaction 
criteria derived above. First, the associated stellar vol- 
ume density estimate is determined for each star within 
the cluster region boundary using the nearest measure- 
ment in our density maps. Note that these are near- 
est neighbor derived volume density estimates, thus they 
may be overestimated. We count the number of stars in 
the cluster regions with densities greater than 10 4 pc~ 3 . 
The mean number of field stars falling within the high 
density areas of the cluster region is estimated directly 
from the Nqm field star models. See Tabled for results 
for each cluster. 

In the two highly embedded clusters presented, 
GGD 12-15 and IRAS 20050+2720, 72% and 91% of the 
inferred cluster region members respectively are in loca- 
tions with stellar densities above 10 4 pc -3 . This result 
suggests that there is ample opportunity for protostellar 
envelope and large T Tauri disk interactions to occur in 
these two clusters even over relatively short timescales 
(< 10 5 yr). Furthermore, while classical T Tauri disk 
scale interactions are less common, they still may occur 
in as many as 10% of the stars in these two clusters over 
the course of ~1 Myr. Indeed, 1 Myr is probably an 
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upper limit on the high stellar density phase lifetime of 
these clusters, as significant gas expulsion is already ap- 
parent. Rega rdless, this estimate is consistent with the 
suggestion of lArmitaee. Clarke, fc Pallal l)2003[) that up 
to 30% of T Tauri disks must be dispersed within 1 Myr, 
possibly due to a combination of close encounters with 
other stars and more gradual evolution and accretion in 
the smallest disks. In contrast, less than 24% of the stars 
in the cluster core boundary of NGC 7129 are at densi- 
ties greater than 10 4 pc -3 . Clearly, there is little chance 
for further disk-affecting interactions i n this cluster, al- 
thoug h this is not entirely unexpected. iGutermuth et all 
(2004) report the disk fraction in the cluster region is 
54%, suggesting that the population is already somewhat 
evolved. Furthermore, spectroscopic evidence suggests 
that the low mass stars in NG C 7129 are between 1.5- 
2.0 Myr old l)HillenbrandllT995|) . 

If we assume cluster asymmetry and high average ex- 
tinction are signs of extreme youth (< 1 Myr) in the 
young clusters presented, then the high density environ- 
ments we observe in GGD 12-15 and IRAS 20050+2720 
provide evidence supporting the idea that protostellar 
envelope and large (10 3 AU) T Tauri disks are likely to 
be affected by close encounters with other members in 
the youngest clusters, and this may lead to significant 
tidal stripping or conversely, accelerated accretion, early 
in their evolution. If the molecular gas mass of the natal 
cloud is indeed adequate to maintain the high density 
configurations we have observed over a 1 Myr timescale, 
a small but measurable fraction of the classical T Tauri 
disks have the chance to be significantly truncated or dis- 
rupted entirely. Once gas expulsion and subsequent dy- 
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namical expansion occurs, as we argue has happened in 
NGC 7129, stellar densities decrease to the level at which 
the disks are allowed to evolve and eventually disperse by 
means of more gradual activity, su ch as photoevaporation 
l)Armitage. Clarke, fc P alla 2003), residual accretion, or 
planet formation. Clearly, determining the timescale of 
bulk gas expulsion is of fundamental importance to un- 
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encounters may have for disk fraction studies, planet 
formation frequency and timescale, and the stellar mass 
spectrum. 
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TABLE 1 

Cluster characteristics compiled from the literature. 



GGD 12-15 IRAS 20050+2720 NGC 7129 



Cluster Distance (pc): 830 a 700 b 1000 c 

Molecular Cloud Mass d (M ): 745 275 400 

Velocity Dispersion d (km s -1 ): 0.92 c 1.07 c 0.62 c 

Far IR Luminosity d (L ): 5680 227 1680 f 



3 Hcrbst & Kacm 
^ Wilking eFai] ti989l) 
IRacinei 11*963? 
IRidee et all 12003*) 

e Con verted from FWHM values (cr = FWHM/2.35) reported in IRidge et~aT1 
<2003T) . 

'Adjusted to account for different adopted distance. 



TABLE 2 

Point Source Detection Completeness Limits 
(90%) 





GGD 12-15 


IRAS 20050+2720 


NGC 7129 


K: 


18.25 


17.20 a 


17.90 a 


H: 


18.85 


18.20 


18.70 


J: 


20.00 


18.70 


19.25 



a K s 



TABLE 3 
Cluster core characteristics. 





GGD 12-15 


IRAS 20050+2720 


NGC 7129 


Cluster-Dominated Region Radius 3, (pc) 


0.29 


0.29 


0.59 


Cluster Core Radius 3, (pc) 


0.24 


0.24 


0.65 


No. of Cluster Region Members 


98 + 10 


91 ± 10 


122 ± 16 


No. of Cluster Region Field Stars 


15 + 2 


58 + 4 


114 ± 8 


Azimuthal Asymmetry Parameter 


1.588 


1.421 


1.207 


Peak Surface Density b (pc — 2 ) 


5910 


6320 


2750 


Mean Surface Density c (pc -2 ) 


371 


344 


91.9 


Peak Volume Density* 3 (pc -3 ) 


2.7 x 10 5 


3.0 x 10 5 


8.6 x 10 4 


Mean Volume Density (pc — 3 ) 


959 


891 


106 



a The cluster-dominated region radius, often more simply called "cluster radius" in the 
text, is used throughout this work in lieu of the cluster core radius. 

b Peak densities are derived from nearest neighbor distances (see Section l3.5l . 

c Mean densities are derived from estimated cluster memberships and cluster-dominated 
region radii. 
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Fig. 1. — Maps of Ak as a function of position over the field of view of our near-IR observations of GGD 12-15. The top map is derived 
purely from the H — K colors of all the stars detected in both bands. The bottom map is the same HK-der'vjed map, updated with values 
derived from the SCUBA 850 fim maps in places where there is a clear underestimate in the H K-derived values. Note that these maps are 
log-scaled and inverted, with locations having Ak < 0.2 in white and Ak > 3.0 in black. The contour levels are Ak = 0.6, 1.2, 2.4, 4.8, 9.6. 
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Fig. 3. — Maps of Ak as a function of position over the field of view of our near-IR observations of NGC 7129. For a detailed description, 
see Fig.0 
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Fig. 4. — i-f-band apparent magnitude histograms for GGD 12-15. The top plot consists of all stars detected within the cluster region 
boundary. The middle plot shows three field star KMH models. The dot-dashed line is the measured off-field KMH, shifted by SAk and 
scaled by the ratio of the areas of th e cluster core and off-fields. The smooth, solid line overplotted on this and the first plot is the field 
star model of Wainscoat ct al. (1993), shifted by the average Ax of the cluster region, for comparison. The solid histogram is the mean 
KMH of the stars within the cluster core boundary from the Ncm = 1000 Monte Carlo field star models generated in Section ^3.21 The 
final plot is the difference of the raw cluster core histogram and a hybrid of the two field star KMHs derived from our observations. The 
mean Monte Carlo model KMH is used up to the limiting magnitude denoted by vertical dashed lines in all three plots. Beyond this limit, 
the scaled and shifted off-field KMH is used. 
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GGD 12-15: Radial Profile 
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Fig. 5. — Azimuthally averaged stellar radial density profiles for the GGD 12-15 region. The solid line histogram on the left shows the 
raw measured radial profile. The overplotted dot-dashed histogram is the mean radial profile from the Nqm = 1000 Monte Carlo field 
star models. The solid line histogram on the right is the contamination-subtracted radial profile. Overplotted on both are smooth solid 
lines showing the exponential plus constant fits, vertical dashed lines marking the radii where the source counts are 3cr ab ove t he field star 
counts, and vertical dotted lines marking the cluster core radii, three times the e-folding lengths from the fits (see Section lij. 31 . 
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Fig. 6. — TOP: Azimuthal distribution histograms of all stars within the cluster core boundary of GGD 12-15. The horizontal dashed line 
marks the mean number of stars per bin, corresponding to an ideally uniform stellar distribution for reference. The dot-dashed histogram 
plotted on the left is the mean azimuthal distribution of the Ncm = 1000 Monte Car lo field star models. BOTTOM: Probability 
distribution histograms for the AAP derived from 10 6 model cluster fields (see Section l3.4l . The left plot shows the raw AAP distribution, 
and the right shows the AAP after subtraction of the recentered mean azimuthal distribution of the field star models. Overplotted on 
both is a dashed line corresponding to a Gaussian of mean |Aj4P 3 i m | and standard deviation <r to demonstrate that the distributions are 
effectively approximated this way. Vertical dotted lines mark the location of the appropriate measured AAP. 
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Fig. 7. — Comparison of three different stellar surface density mapping methods for the three clusters presented: G GD 12-15, 
IRAS 20050+2720, and NGC 7129. In the left column are the N = 5 stellar density maps used throughout this work (see Section |3~51 . The 
maps in the central column are made with the same method, but the densities at each gri d point are calculated using the distance to the 
nearest N = 10 stars. Maps derived using Gaussian kernel smoothing (Gomez et al. 1993) with a smoothing width of 8" are in the right 
column. All maps are 5' on a side. 
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Fig. 8. — Plot of all i-f-band detections brighter than K = 18 in GGD 12-15. The box marks the off-field u sed for characterization of 
local field star contamination. The circle marks the cluster region boundary as determined in Fi g. [51 (see Section 13.31 . The overlapping 45° 
position angle bins used to make the Nyquist sampled azimuthal distribution histograms in Fig.jfcjare marked by the alternating black and 
gray wedge outlines. 



22 



Gutermuth et al 




RA (J20D0) 

Fig. 9. — JHK log-scaled color image of the GGD 12-15 region from FLAMINGOS on MMT overlaid with 850 /im dust emission 
contours from SCUBA on JCMT. The contour levels start at 1.75 mjy pixel -1 (2" X 2" pixels), with successive contour levels at double 
the previous level. The blue circle shows the cluster region boundary determined in Fig. [F] 
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RA (J20DO) 

Fig. 10. — 850 fim SCUBA map of GGD 12-15 in negative grayscale, scaled to show faint filamentary structure. Stellar surface density 
contours measured from all if -band point sources brighter than K = 18 are overlaid. Contours begin at 435 pc~ 2 (5a above median field 
star density), and increase at an interval of 750 pc -2 . The vectors denote the location and extent of the stellar surface density map cuts 
plotted in Fig. 1111 Each vector is 1 pc in length at the assumed distance of 830 pc, approximately the distance a star moving lfcms -1 will 
travel in 1 Myr. 
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Fig. 11. — The black lines in the above plots are stellar surface density map cuts for GGD 12-15. The vectors plotted in Fig. UTil mark 
the position, orientation, and direction of increasing offset from cluster center. Note the strong central peak, multiple sub-peaks, and lack 
of uniform circular sym metry . The gray line marks the 5cr above median field star density level for comparison, which corresponds to the 
first contour level in Fig. I1UI 



12G 
Si BD 

i 4 



Spatial Distribution Analysis in Young Clusters 

P£5 ^QdKi^Sff XHH; Fn* Clmriar Pugrnj 14B ttcs 



25 





I 






- 


! 








I 








i 


J 






1 
l 

1 , 

; — — ~i 


J 

/ 


-i 



iz ■ n in 



13- 



30 



IIP 
10 D 
HP 



2D 
P 



Wfl3: A, - 1.90 

K 

FLAW: 4\ - 1.2S 
p- 



Ita 




ID 12 14 1fl 14 

AffMPHft MflptUrif 

IRAS 2JMDO+ZT5D KMH: CM-nubU ClLritnr: Bl 




H lfl 

Fig. 12. — iC s — band apparent magnitude histograms for IRAS 20050+2720. For a detailed description, see Fig. ITI 
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IRAS 20050 + 2720: Radial Profile CM-sub'd Radial Profile 
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Fig. 13. — Azimuthally averaged stellar radial density profiles for the IRAS 20050+2720 region. For a detailed description, see Fig. 151 



Spatial Distribution Analysis in Young Clusters 
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Fig. 14. — Azimuthal distribution histograms and AAP probability distribution histograms for IRAS 20050+2720. For a detailed 
description, see Fig. 151 
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20050+2720: all aaurcBa with K. z 1 B.O 




Fig. 15. — Plot of all if s -band detections brighter than K s = 16 in IRAS 20050+2720. The box mar ks th e off-field us ed fo r characterization 
of local field star contamination. The circle marks the cluster region boundary as determined in Fig. 1131 (see Section 13.31 . The overlapping 
45° positional angle bins used to make the Nyquist sampled azimuthal distribution histograms in Fig. 1141 are marked by the alternating 
black and gray wedge outlines. 
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Fig. 16. — JHK S log-scaled color image of the IRAS 20050+2720 region from FLAMINGOS on MMT overlaid with 850 /im dust emission 
contours from SCUBA on JCMT. The contour levels start at 1.75 mjy pixel -1 (2" X 2" pixe ls), with successive contour levels at two times 
the previous level. The blue circle shows the cluster region boundary determined in Fig. 1151 
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Fig. 17. — 850 fim SCUBA map of IRAS 20050+2720 in negative grayscale, scaled to show faint filamentary structure. Stellar surface 
density contours measured from all i\s-band point sources brighter than K a = 16 are overlaid. Contours begin at 1450 pc — 2 (5cr above 
median field star density), and increase at an interval of 750 pc — 2 . The vectors denote the location and extent of the stellar density map 
cuts plotted in Fig. 1181 Each vector is 1 pc in length at the assumed distance of 700 pc, approximately the distance a star moving lfcms -1 
will travel in 1 M yr. The letters A, B, and C denote the contour groups that correspond to the three subclusters given those names by 
IChen et alj Il997t) . 
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Fig. 18. — The black lines in the above plots are stellar surface density map cuts for IRAS 20050+2720. The vectors plotted in Fig. El 
mark the position, orientation, and direction of increasing offset from cluster center. Note the multiple strong peaks and deviation from 
circular sy mmet ry. The gray line marks the 5<r above median field star density level for comparison, which corresponds to the first contour 
level in Fig. El 
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Fig. 19. — K a — band apparent magnitude histograms for NGC 7129. For a detailed description, see Fig. El 



Spatial Distribution Analysis in Young Clusters 



33 
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Fig. 20. — Azimuthally averaged stellar radial density profiles for the NGC 7129 region. For a detailed description, see Fig. 151 
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Fig. 21. — Azimuthal distribution histograms and AAP probability distribution histograms for NGC 7129. For a detailed description, 
see Fig. ITT1 
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NQC 7129: all aourcaa with K a < 17.5 
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Fig. 22. — Plot of all X s -band detections brighter than K B = 17.5 in NGC 7129. The box marks the off-field us ed fo r characterization 
of local field star contamination. The circle marks the cluster region boundary as determined in Fig. 1201 (see Section The overlapping 
45° positional angle bins used to make the Nyquist sampled azimuthal distribution histograms in Fig. 1211 are marked by the alternating 
black and gray wedge outlines. 
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Fig. 23. — JHK 3 log-scaled color image of the NGC 7129 region from FLAMINGOS on MMT overlaid with 850 /an dust emission 
contours from SCUBA on JCMT. The contour levels start at 1.75 mjy pixel -1 (2" X 2 " pi xels), with successive contour levels at double 
the previous level. The blue circle shows the cluster region boundary determined in Fig. 121 II 
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Fig. 24. — 850 /im SCUBA map of NGC 7129 in negative grayscale, scaled to show faint filamentary structure. Stellar surface density 
contours measured from all A'a-band point sources brighter than K s = 17.5 are overlaid. Contours begin at 475 pc -2 (5cr above median 
field star density), and increase at an interval of 750 pc -2 . The vectors denote the location and extent of the stellar density map cuts 
plotted in Fig. 1251 Each vector is 1 pc in length at the assumed distance of 1 kpc, approximately the distance a star moving ltms' 1 will 
travel in 1 Myr. 




Fig. 25. — The black lines in the above plots are stellar surface density map cuts for NGC 7129. The vectors plotted in Fig. l24l mark the 
position, orientation, and direction of increasing offset from cluster center. Note the lack of a cohesive peak and overall diffuse stru ctur e. 
The gray line marks the 5cr above median field star density level for comparison, which corresponds to the first contour level in Fig. 1241 
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Fig. 26. — Estimates for collisional timescales for protostellar envelopes (R = 10 3 AU, solid line) and T Tauri disks (R = 10 2 AU, dashed 
line) as a function of stellar volume density. 
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TABLE 4 

Nearest Neighbor. Volume Density Environments by Cluster 



GGD 12-15 IRAS 20050+2720 NGC 7129 



No. of Stars at > 10 4 pc" 3 : 
Est. No. of Field Stars: 
No. of Members at > 10 4 pc~ 3 : 
Fraction of total in core: 



74 100 35 
3 17 6 

71 83 29 

72.4% 91.2% 23.8% 



